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The general solution of the gravitational field equations for a full causal bulk viscous stiff cosmo- 
logical fluid, with bulk viscosity coefficient proportional to the energy density to the power 1/4, is ob- 
tained in the fiat Friedmann-Robertson- Walker geometry. The solution describes a non-infiationary 
Universe, which starts its evolution from a singular state. The time variation of the scale factor, 
deceleration parameter, viscous pressure, viscous pressure-thermodynamic pressure ratio, comoving 
entropy and Ricci and Kretschmann invariants is considered in detail. 
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O \ Dissipative bulk viscous type thermodynamical processes are supposed to play a crucial role in the dynamics 
and evolution of the early Universe. There are many processes capable of producing bulk viscous stresses in the 
early cosmological fluid, like interaction between matter and radiation, quark and gluon plasma viscosity, different 
components of dark matter, particle creation, strings and topological defects etc. Traditionally, for the description 
of these phenomena the theories of Eckart and Landau and Lifshitz were used. But Hiscock and Lindblom Q 
have shown that the Eckart-type theories suffer from serious drawbacks concerning causality and stability. Regardless 
^ of the choice of equation of state, all equilibrium states in these theories are unstable and in addition signals may be 
propagated through the fluid at velocities exceeding the speed of light. These problems arise due to the first order 
^n nature of the theory, i.e., it considers only first-order deviations from the equilibrium. The neglected second-order 
^— ^ terms are necessary to prevent non-causal and unstable behavior. 

I , A relativistic consistent second-order theory was found by Israel Q and developed into what is called transient 
■ or extended irreversible thermodynamics |^, 0], 0- Hiscock and Lindblom have shown that these second order 
. theories are free of the pathologies of the Eckart-type theories. Therefore, the best currently available theory for 
' analyzing dissipative processes in the Universe is the full Israel-Stewart-Hiscock causal thermodynamics. 

Exact general solutions of the field equations for flat homogeneous Universes filled with a full causal viscous fluid 
I ■ source for a power-law dependence of the bulk viscosity coefficient on the energy density have been obtained in 0- 
^jrjl p^. The evolution of a homogeneous and isotropic dissipative fluid in the truncated Israel-Stewart theory has been 
• • . analyzed, by using dynamical systems methods, by Di Prisco, Hcrrera and Ibanez [Tll |. They have found that almost 
. ^ ' solutions inflate, and only a few of them can be considered physical, since they do not satisfy the dominant energy 

, condition. 

^ It is the purpose of the present paper to obtain the general solution of the gravitational field equations for the 

case of a bulk viscous cosmological fluid in a flat Friedmann-Robertson- Walker type geometry, obeying the Zeldovich 
(stiff) equation of state, with the bulk viscosity coefficient ^ proportional to the energy density p to the power 1/4, 
^ cx p^/^. In this case the solution of the gravitational field equations can be represented in an exact parametric 
form, in terms of Bessel functions. The behavior of the scale factor, deceleration parameter, viscous pressure, viscous 
pressure-thermodynamic pressure ratio, comoving entropy and Ricci and Kretschmann invariants is also considered. 
The solution describes a non-inflationary Universe, which starts its evolution from a singular state. 

The present paper is organized as follows. The basic equations describing the dynamics of the viscous Friedmann- 
Robertson- Walker (FRW) Universe are written down in Section II. The general solution of the field equations is 
obtained in Section III. In Section IV we discuss and conclude our results. 
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II. GEOMETRY, FIELD EQUATIONS AND CONSEQUENCES 

The energy-momentum tensor for a bulk viscous cosmological fluid is given by 

Tt = {p + P + Tl)u,u^-{p + Il)St (1) 

where p is the energy density, p the thermodynamic pressure, 11 the bulk viscous pressure and Ui the four velocity 
satisfying the condition u^m* = 1. We use units so that SttG — c — 1. 

For a flat Friedmann- Robertson- Walker Universe with a line element given by ds^ — dt^ — a'^{t) [dx^ + dy^ + dz^), 
filled with a bulk viscous cosmological fluid, the gravitational field equations and the continuity equation T*^ — 
(with ; denoting the covariant derivative with respect to the metric) imply 

3ij2 = p,2H + 3H^ = -p-Il,p + 3{p + p)H = -SHU, (2) 

where H = d/a is the Hubble parameter. The causal evolution equation for the bulk viscous pressure is given by [l^ 

Ttl + n^-3^H-Kll(^3H+^-i-^y (3) 

where T is the temperature, ^ the bulk viscosity coefficient and t the relaxation time. Eq. Q arises as the simplest 
way (linear in H) to satisfy the H theorem ( i.e., for the entropy production to be non- negative, S\ = iP /£^T > 0, 

where 5* — eN^ ~ ^j^u'' is the entropy flow vector, e is the entropy per particle and — nu^ is the particle flow 
vector) 0, 0. 

In order to close the system of equations Q we have to give the equation of state for p and specify T, ^ and r. As 
usual, we assume the following phenomenological laws T2]: 

P={l-l)p,i = ap',T = /3p^ r - ^p-' = ap'-\ (4) 

where 1 < 7 < 2, and a>0, /3>0, r>0 and s > are constants. Eqs. I@J are standard in cosmological models 
whereas the equation for t is a simple procedure to ensure that the speed of viscous pulses does not exceed the speed 
of hght. 

The requirement that the entropy is a state function imposes in the present model the constraint J' = (7 — 1) /7 Q 
so that < r < 1/2 for 1 < 7 < 2. 

The growth of the total comoving entropy E over a proper time interval {to,t) is given by |12| |: 

E(t) - E (to) = -Sfc^^ [ UHa^T^^dt, (5) 

J to 

where ks is the Boltzmann constant. The Israel-Stewart-Hiscock theory is derived under the assumption that the 
thermodynamical state of the fluid is close to equilibrium, that is the non-equilibrium bulk viscous pressure should 
be small when compared to the local equilibrium pressure |n| << p = (7 — 1) p 13]. If this condition is violated then 
one is effectively assuming that the linear theory holds also in the nonlinear regime far from equilibrium. For a fluid 
description of the matter, the condition ought to be satisfied. 

To see if a cosmological model inflates or not it is convenient to introduce the deceleration parameter q = 
d {H'^) /dt — I ^ [p + Zp + 3n) /2p. The positive sign of the deceleration parameter corresponds to standard decel- 
erating models, whereas the negative sign indicates inflation. 

With these assumptions the evolution equation for the fiat homogeneous causal bulk viscous cosmological models 

is 

ii + 2,HH + •i^-'a-^H'^-^'H - (1 + r)H-^H^ + ^ (^ - 2) + ^7^4-2. = q. (6) 

4 2a 

By introducing a set of non-dimensional variables h and 6 by means of the transformations H — a^h^t — 3^^, 

s 1 

with ao — (■^^) I s 7^ 51 E^nd using the expression of r as a function of 7, Eq. ^ takes the form 

d^h r on ,^l dh , , , f dh\^ 2r — 1 „ 1 o/o 

^ + + '^'''-''] ^ - (1 + f^-' (^) + Y^^' + T-r' ' ^ ' = 

By denoting n — {1 — 2s) / (1 — r) and changing the variables according to h — y^/'-^^^\ri — J y^^^^^^^dd, Eq. 
becomes: 



^ + (2 + y")^ + (2r-l)j/ + y"+i=0. (8) 
ayy-^ df] 
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III. GENERAL SOLUTION FOR A STIFF COSMOLOGICAL FLUID WITH s = 1/4 

The stiff cosmological fluid, with equilibrium pressure equal to the energy density, is supposed to give an accurate 
description of the very early phases of the evolution of the Universe. It corresponds to the values of the parameters 
7 = 2 and, consequently, r = 1/2. We assume that the bulk viscosity coefficient ^ is of the form ^ = ap^/**, with 
s = 1/4. Hence it follows that n = 1. Therefore the nonlinear second order ordinary differential equation ((SJ becomes: 



ar;^ ar] 



(9) 



A particular solution of Eq. for n = 1, s = 1/4 has been obtained in |1^. The general solution of Eq. ^ is 
given by 



2Cie-^ [Ji {Cie-^) + C2Y1 (Cie-")] 
C2I0 {Cie-v) + Jo (Cie-") 



(10) 



where Ci and C2 are arbitrary integration constants, Jk (x) and Yk {x) are the Bessel functions of the first and second 
kinds, respectively. They satisfy the Bessel differential equation 



(11) 



Hence the general solution of the gravitational field equations for a Zeldovich causal bulk viscous fluid filled fiat 
FRW Universe, with bulk viscosity coefficient ^ = cep^^'^, can be represented in the following exact parametric form, 
with 7] taken as parameter: 



t (??) - to 



3ao J l^Ir]) ' ^ ^ ' ^ ^° i^v] .Pil) ^ P iv) = Sagy^ (77) , 



(12) 



3 dy 
y iv) drj 



-1,H(77) = -6a^y^ (7?) 



drj 



= 2 



dlny 



drj 



(13) 



S (77) = S (r/o) + 



4v^, 



aoal 



kBl3 



y{v) 



( \ dy 



exp (2?7) dr], 



(14) 



where Oq, to and S (770) are constants of integration and ao — 2 /3a. 

The singular or non-singular character of the solution for all times t > can be checked from the finite (infinite) 
character of the Ricci invariant RijR^^ and Kretschmann scalar RijkiR'-'^^ , given by 



I = RijR'^ = 12 



a I a 
a \a 



3(^) +3^+,^ 
077 / a?7 



(15) 



J — RijkiR^'' 



kl 



12 



12a4y6 



drj J drj 



(16) 



IV. DISCUSSIONS AND FINAL REMARKS 



The evolution of the causal bulk viscous Zeldovich fiuid filled fiat Universe starts generally its evolution from a 
singular state, as can be seen from the singular behavior of the invariant RijR^^ , represented, for different values of 
the integration constants Ci and C2, in Fig.l. 



Generally, the dynamics of the Universe depends on the numerical values of the constants Ci and C2 . The behavior 
of the scale factor is presented in Fig. 2, for some specific values of the integration constants. The evolution is 
expansionary for all times. 
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FIG. 1: Time evolution of the invariant RijR^-' for different values of the integration constants Ci and C2: Ci = 1 and C2 = — 1 
(solid curve), Ci = 1.1 and C2 = —0.9 (dotted curve) and Ci = 1.1 and C2 = —1 (dashed curve). 
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FIG. 2: Time variation of the scale factor a for different values of the integration constants Ci and C2: Ci = 1 and C2 = — 1 
(solid curve), Ci = 1.1 and C2 ~ —0.9 (dotted curve) and Ci = 1.1 and C2 ~ —1 (dashed curve). 



The variation of the energy density of the cosmological fluid is represented in Fig. 3. 
At the initial moment the scale factor is zero, while the energy density tends to infinity. 

The dynamics of the deceleration parameter, shown in Fig. 4, indicates, for the chosen range of the integration 
constant, a non-inflationary behavior for all times, with q > 0,yt > 0. 

The bulk viscous pressure 11 is negative during the cosmological evolution, 11 < 0,Vt > 0, as expected from a 
thermodynamic point of view. In the large time limit the viscous pressure tends to zero. In the same limit the bulk 
viscosity coefficient also becomes negligible small. During the cosmological evolution a large amount of comoving 
entropy is produced, with the entropy E increasing in time and tending in the large time limit to a constant value. 
The time variation of the entropy is presented in Fig. 5. 

The ratio of the bulk and thermodynamic pressures, |n/p| is presented in Fig. 6. 

For the chosen values of the integration constants the condition Il/p does not hold for all times, and therefore 
the model is thermodynamically consistent only for a finite period of time. Therefore the solution represented by 
Eqs. H12|l - H14() can consistently describe only a well-determined, non-inflationary period of the early dynamics of the 
super-dense, matter dominated post-inflationary era, when, as expected, the bulk viscous dissipative effects play an 
important role. 
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FIG. 3: Evolution of the energy density p for different values of the integration constants Ci and C2: Ci = 1 and C2 
(solid curve), Ci = 1.2 and C2 ~ —0.8 (dotted curve) and Ci = 1.1 and C2 ~ —1.8 (dashed curve). 
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FIG. 4: Dynamics of the deceleration parameter q for different values of the integration constants Ci and C2: Ci = 1 and 
C2 = —1 (solid curve), Ci = 1.2 and C2 = —0.8 (dotted curve) and Ci = 1.1 and C2 = —1.8 (dashed curve). 
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FIG. 5: Behavior of the comoving entropy E for different values of the integration constants Ci and C2: Ci = 1 and C2 = — 1 
(solid curve), Ci = 1.05 and C2 = —0.9 (dotted curve) and Ci = 1.08 and C2 = —0.8 (dashed curve). The numerical parameters 
have been normalized so that 4%/3QoOo/fcs/3 = 1. 
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FIG. 6: Variation with time of the ratio S for different values of the integration constants Ci and C2: Ci = 1 and C2 = — 1 
(solid curve), Ci = 1.1 and C2 = —0.9 (dotted curve) and Ci = 1.2 and C2 = —0.9 (dashed curve). 



